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Abstract 


The so-called equation of motion method is useful to obtain the explicit 
form of the eigenvectors and eigenvalues of certain non self-adjoint bosonic 
Hamiltonians with real eigenvalues. These operators can be diagonalized 
when they are expressed in terms of pseudo-bosons, which do not behave 
as ordinary bosons under the adjoint transformation, but obey the Weil- 


Heisenberg commutation relations. 

1 Introduction and preliminaries 

In conventional formulations of non-relativistic quantum mechanics, the Hamil¬ 
tonian operator is self-adjoint. However, certain relativistic extensions of quan¬ 
tum mechanics lead to the consideration of non self-adjoint Hamiltonian operators 
with a real spectrum. This motivated an intense research activity, both on the 
physical and mathematical level (see, e.g. [2l IH El El [71 [HI IS] and their references). 

Throughout, we shall use synonymously the terms Hermitian and self-adjoint. 
Denote by the Hilbert space of square integrable functions in two real 

variables, endowed with the standard inner product 
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For 


V = |/(x,|/)e ; f(x,y) is a polynomial in x,yj (1) 

which is a dense domain in [3], let a, fe : "D —)■ P be bosonic operators 

defined, respectively, by 


Id Id 

“ = ^+ 291 - '’ = "+ 2 %' 


Consider also their adjoints a*,b* ■. V ^ V 


* ^ ^ U* ^ ^ 

a =x -—, b =y -—. 

2dx 2dy 


We recall that, conventionally, a, b are said to be annihilation operators, while 
a*,b* are creation operators. It is worth noticing that these operators are un¬ 
bounded. Moreover, P is stable under the action of a, b and of their adjoints, and 
they satisfy the commutation rules (CR’s), 


[a,a*] = [b,b*] = l, 


( 2 ) 


where 1 is the identity operator in "D. (This means that aa*f — a*af = bb*f — 
b*bf = f for any / G P.) Furthermore, 


[a, b*] = [b, a*] = [a*, b*] = [a, b] = 0. (3) 

As it is well-known, the canonical commutation relations ([2]) and ([2]) characterize 
an algebra of Weil-Heisenberg (W-H). Moreover, the following holds. 


a<ho = ^*^0 = 0, 


for $0 = e ^ -p, a so-called vacuum state. The set of vectors 


= m,n>0}, (4) 

constitutes a basis of "H, that is, every vector in can be uniquely expressed 

in terms of this vector system, which is complete, since 0 is the only vector 
orthogonal to all its elements. 

The main goal of this note is to investigate spectral properties of certain non 
self-adjoint operators which are expressed as quadratic combinations of bosonic 
operators. 


2 Non self-adjoint Hamiltonian constructed in 
terms of sr(1, 1) generators 

A linear operator acting on a finite dimensional Hilbert space which is non 
self-adjoint and has distinct real eigenvalues, is similar to its adjoint operator. 
Concretely, the following holds, having in mind that the spectrum of a finite 
matrix is the set of its eigenvalues. 
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Theorem 2.1 Let H be an n x n complex non self-adjoint matrix with distinct 
real eigenvalues. Then H and H* have a common spectrum, a complete systems 
of eigenvectors and they are similar, that is, there exists a unitary matrix S such 
that H* = SHS~^. Moreover, if Heft = X^cf, H*fji = Xifji, i = 1,... ,n, then the 
eigenvectors may be normalized so that = 5ij. 

Proof. Since the eigenvalues of H are real, they coincide with those of H* . Let 
cr(iL) = cr{H*) = {Ai,..., A„}, and let <fi, iji, i = 1,... ,n be such that 

H(j)i = Xicfi, H*i)i = Xi'ijji. 

Further, as the eigenvalues are mutually distinct, the corresponding eigenvectors 
are linearly independent, and so the eigensystems {'iji}, {(fj} corresponding to 
the eigenvalues A,, i = 1,2,... ,n, are complete. 

Let S be dehned by xji = S(j)i. Thus, S is unitary and 

H*S(t)i = XiSefi = SHcfi, i = l,...,n, 

implying that H and H* are unitarily similar. 

Having in mind that 

the orthonormality relation follows. ■ 

The following question naturally arises. Do the above properties survive in 
the inhnite dimensional setting? We answer this question for the model described 
by the linear operator H : —>■ dehned as 

H = a*a + bb* + /3{a*a — bb*) + 'y{a*b* — ah), /?,7 G M, (5) 

which treats a system of two interacting bosons. It is obvious that H is non 
self-adjoint. 

Consider the class C of unbound linear operators H acting on the inhnite 
dimensional Hilbert space TL and satisfying the following properties: 

(I) H is non self-adjoint, 

(H) H has real eigenvalues, 

(HI) H and H* are isospectral. 

(IV) The associated eigenvectors form biorthogonal bases of TL. 

(V) There exists a unitary transformation S on TL such that H* = SHS~^. 

Throughout we prove that iL in ([5]) belongs to the class C, and so, in the 
eigenvectors bases, H and H* have a diagonal representation. 
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3 The equation of motion method for pseudo- 
bosonic operators 

Let us consider a non self-adjoint Hamiltonian H [H ^ H*) which is expressed 
in terms of unbounded bosonic operators Oi,..., a„, a^, ..., a*, that is, linear 
operators acting on V and satisfying the CR’s, 

[ai,a*] = 5ij, [ai,aj] = 0 , [a*, a*] = 0 . ( 6 ) 

As usual, 6ij denotes the Kroenecker symbol (which equals 0 for i 7 ^ j and 1 for 

i = j)- 

Let V be the linear space spanned by oi,..., a„, a^,..., a* and let / G V, 
that is, / is a linear combination of the operators ai,..., a„, aj,..., a*. Having 
in mind that 

and their adjoint relations, it is clear that [if, /] G V. Obviously, the operator 
[if, /] : V —)■ V is a linear operator on / and next we prove that its eigenvalues 
are closely related with those of H . 

Theorem 3.1 Let H he a non self-adjoint operator acting on ji with real eigen¬ 
values expressed in terms of boson operators satisfying Let / G V and A G M 
be such that 

W, f] = A/. (7) 

Suppose A G M and xj ^ TL are, respectively, an eigenvalue of H and an asso¬ 
ciated eigenvector. If fxjj ^ 0, then also A -|- A is an eigenvalue of H and fxp a 
corresponding eigenvector. Moreover, there exists 0 7 ^ 5 ^ G V such that 

[H,g] = -Xg. 

Proof. Under the hypothesis, we have Hip = Axp. Then, 

HfxP = [H, f]xp + fHfj = (A + X)fxp. 

Since A and A -|- A are both real, so is A. Since [H*,f*] = —Xf* and H and 
H* have the same eigenvalues, it follows that there exists 0 7 ^ 5 ^ G V such that 
[H, g] = —Xg. Observe that the operator [H, /] : V —)■ V is represented by some 
2n X 2n matrix T, while the operator [H*,f] : V ^ V is represented by the 
transconjugate matrix T*, and so the eigenvalues of T and T* coincide since they 
are real. ■ 

If A > 0, / is called an excitation operator [deexcitation operator \i A<0). 
The so called equation of motion method (EMM) consists in the solution 
of eq. o. with the determination of the eigenvalues of H and corresponding 
eigenvectors. 
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Let 0 7 ^ G V be such that [H^g] = \'g. The Jacobi identity yields 

l|if,/],9l + [l/,9l,-ffl + [|9,-ffl,/] = 0. 

Now, [[/, = 0, since [/, ( 7 ] is a multiple of the identity. Thus, 

(A + X')[f,g] = 0. 

As a consequence, either A = —A' or [/, ^f] = 0. If [f,g] 7 ^ 0, which happens only 
if A = —A', / and g may be normalized so that [f,g] = 1. Since f* 7 ^ g, the 
operators obtained by this procedure do not describe “ordinary” bosons. In this 
case, it may be seen that there exists 0 7 ^ J/q G U such that either /Tq = 0 
or ( 7 T 0 = 0. Assume that /Tq = 0. We write g = P and it can be checked 
that the vectors system {/I "Tq : n > 0} constitutes a complete basis of 
There also exists Tg G U such that g*'^Q = = 0 and the vector system 

{/* "J/q : n > 0} is an eigenvector system for H*, which constitutes a complete 
basis of The bases 

: n > 0 } and {f* ’"T'g : n > 0 } 

are biorthogonal: 

Thus, the operators /, are said to describe pseudo-bosons for if, while the 
operators p*, f* describe pseudo-bosons for H* (cf. Bagarello [HIS]). If = 0, 
with J/q G V, an analogous discussion takes place. 

4 H belongs to the class C 

Throughout we prove that H in eq. ([5]) belongs to the class C. As already 
observed, H satishes (I). 

(II) H has real eigenvalues. 

The operator H is quadratic in the bosonic operators and we may explicitly 
hnd its eigenvalues by the equations of motion method, as follows. We look for 
a linear combination of the bosonic operators such that its commutator with H 
satishes the proportionality condition, 

[H, {xaa* + Xbb* -h Vaa -h ytb)] = X{xaa* x^b* + i/aa ybb), A G M. 

This leads to the following system of linear equations in Xa, Xb, ya, yt and A: 

(1 /3)Xa - jyb = Xxa, 

(1 - f3)xb - ^ya = Axfe, 

-{1 + (5)ya-lXb = Xya, 

-(1 - I3)yb - 7 X„ = Xyb. 
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Solving this linear system, we readily obtain A = ±/5 ± + 7 ^. The vectors 

{xa, Xb, Va, VbY' and the parameters A are, respectively, the eigenvectors and eigen¬ 
values of the 4x4 real symmetric matrix 



1 + /3 

0 

0 

-7 

0 

1-/9 

-7 

0 

0 

-7 

-1-/? 

0 

-7 

0 

0 

—1 -|- /3 


The eigenvalues of T and associated eigenvectors are given by 


Ai = -/3 - a/1 72, vi = A/'(0, -1 a/1 + 72,7, 0)^, 

A2 =/9 - a/TT^, V 2 = -h v^l + 72,0,0,7)^, 

As = -/^ + Vl + 7^ = ■A/'(0,1-h ^/TT^,-7,0)^, 

A4 =/^ + a/TT^, Va = A^(l a/1 4- 72, 0, 0, . 

The determination of the normalization constant M is postponed and done ac¬ 
cording to future convenience. 

Consider now the operators 

c = A^((-l + \/l + 72)6* 7a), 

d = ^/l -h 7^)a* -F 76), 

= A/’((l -h \/l -h 72)6* — 7a), 

= AA((1 -F a /1 72)0* - 76), 

also dehned on D. Notice that 7 ^ c* and ^ d*. 

Let us take 

M = ^27a/ 1 72 j 

Then, the operators c^, d^, c, d, satisfy the CR’s of a Weil-Heisenberg algebra, 

[c, c^] = [d, d^] = 1 , 

vanishing all the remaining commutators between them. Moreover, H can be 
written as 

H = /9(c^c — d^d) -I- \/l + 72(c^c -|- dd^). 

Next, we show that the Hamiltonian presents a diagonal form when it is expressed 
in terms of the operators c^, d^, c, d. 

Having in mind that 


$0 = exp {-{x‘^ + y^)) e V 
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satisfies a<f)o = ^*^*0 = 0 , it may be easily verified that 

d'o = exp(-aa* 6 *)<ho = exp {x^ + 2 /^) - , 

\ 1 — a"' 1 — J 

where a = 7/(1 + a/I + 7 ^) = (—1 + + 7 ^)/ 7 , satishes 

cd'o = dd'o = 0. 

Of course, \l/o G spanP. Let 

{d'm.n = : m,n>0}. 

From the CR’s of a Weil-Heisenberg algebra, it follows that 

H^m,n = (^a /1 + 7 ^ + + Vl + 7^) + ’^(-/^ + Vl + 7^)) 

Henceforth, the eigenvalues of H are given by 

Em,n = \/l + 7^ + m{l3 + v^l + 7 ^) + n(-/? + V"! + y^). ( 8 ) 

Thus, (II) holds. Moreover, the associated eigenvectors to Emn are 

^m,n = m,n > 0 . 

These eigenvectors constitute a basis of Ti. Further, it can be shown that this 
basis is complete [3]. 

In order to prove that id* and El have the same eigenvalues, let us consider 


the state vector 

Tq = exp(aa*6*)<I)o 

Clearly, 

c^T'o = = 0. 

It follows that the state vectors 


= m,n 


are eigenvectors of H*, 

= (71+7+m(/3+\/r+7)+ n{-i3+ yr+v)) 

Hence, H* and H are isospectral and (HI) holds. 

Next, we show that the basis constituted by the eigenvectors of H is orthog¬ 
onal to the basis of the eigenvectors of H*. The basis 
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is orthogonal to the basis ■ m,n> 0}, as we have 

= m\n\5^p5n,g{^o,%)■ 

B.ence^oTth., (IV) is satished. 

Finally, (V) holds. The nnitary transformation 

S = exp + 6*6) j , 

acts on a and 6 as follows: 

a ^ ia = SaS~^, b ^ ib = SbS~^, 
and so S takes H into H* , that is, 

H* = SHS-\ 


Hence, H satishes (V). 


Some considerations are in order. 

1 . In the above considered eigenvectors system, the matrix of H has the 
following representation 


p 


0 

0 

0 



0 

(3 + 2 p 

0 

0 



0 


0 

2 ( 3 + 3 p 

0 


®I 

0 


0 

0 

3(3 + 4 p 



*- 


0 

0 

0 


0 



0 - 

■(3 + p 

0 


0 

+/ 

(8) 

0 

0 2(- 

-(3 + p) 


0 



0 

0 

0 

3 {—( 3 -\rp) ... 


where / is the semi-inhnite identity matrix, p = + 7^ and ® denotes the usual 

Kroenecker product. Hence, the matrix of H is the following diagonal matrix 


00 


© 


k {(3 + p)+ p 
0 
0 


0 0 
kiyjd -\- p) — /3 + 2p 0 

0 /c (/3 + p) — ‘ 2/3 + 3 p 


2 . In the consstructed eigenvectors bases, H and H* have a diagonal repre¬ 
sentation. 











3 . The integers m,n are the number of pseudo-bosons of type c^,d^, re¬ 
spectively. Since, however, 7^ c*, d} ^ d*, they do not describe “ordinary” 
dynamical bosons and they are said to describe dynamical pseudo-bosons BE]. 

Remark that, although H is non self-adjoint, it has real eigenvalues and a 
complete system of eigenvectors for any 7, /3 G M. If we replace 7 G M by 
i\, H becomes self-adjoint, but then it only has real eigenvalues and a com¬ 
plete system of eigenvectors if I7I < 1 . The existence of complex eigenvalues for 
I7I > 1 indicates that the system under consideration is not stable, or that H is 
not bounded from below, violating conditions which, in general, are required on 
physical grounds. 

5 Invariant subspaces 

In this Section, we provide an alternative approach to determine the eigen¬ 
values and eigenvectors of H. 

We notice that the operator H is expressed as a linear combination of the 
generators of the sm(1, 1) algebra 

a*a + bb*, ab, a*b*, 

and of the operator 

a* a — b*b, 

called, for simplicity, Casimir operator, since it is closely related to the actual 
Casimir operator which is given by 

C := {a*a + bb*f - 2 a*b*ab - 2 aba*b* = {a*a - b*b - l){a*a - b*b + 1). 

This is the key fact for the alternative procedure for the spectral analysis in this 
Section. 

Since any / G may be expanded in the basis {<hmn : m,n > 0 }, we 

may identify with Ti = span{<hmn ■ m,n> 0 }. It is convenient to express 

T-L as a. direct sum of eigenspaces of the Casimir operator, 

+ CXD 

'H= ^ T-l-k, "Hfc = span{<hm„;m - n = fc,m,n > 0}. 

k=—oo 

Notice that eigenspaces of the Casimir operator are invariant subspaces of the 
generators of the sm(1, 1 ) algebra, and consequently of H. In Tik, the generator 
a*b* is matricially represented by 


r 0 

0 

0 

0 ...1 

\/\k\ 1 

0 

0 

0 ... 

0 

72(|A^|+2) 

0 

0 ... 

0 

0 

V'3( A; 3) 

0 ... 
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while the generator ab is represented by 


A_ 





and the generator a*a + bb* is represented by 


Aq — 


\k\ + l 0 0 0 

0 |fc|+3 0 0 

0 0 |fc|+5 0 

0 0 0 \k\ + 7 


The following CR’s of the sn(l, 1) algebra are satisfied 

[A-^Aj^\ = Aq, [^40,^+] = 2R+, [Aq, A_] = —2A_. (9) 

The operator A^ (^-) is said to be a raising (lowering) operator. That is, if is 
an eigenvector of Aq, so that = A$, then A+d) is an eigenvector associated 
with an upward shifted eigenvalue, 


AoA+^ = (A + 2)A+<f). 

Similarly, A_<h 7^ 0 is an eigenvector associated with a downward shifted eigen¬ 
value, 

AoA_<h = (A - 2)A_<h. 

The spectrum of Aq is bounded from below and the eigenvector $0 = (1, 0, 0,.. .)^ 
satisfies A_<f)o = 0, being called a lowest weight state. A set of eigenvectors associ¬ 
ated with eigenvalues of Aq, which are positive, is obtained by acting successively 
with A_|_ on $0. We observe that Aq = Aq and A+ = A*_. It is interesting 
that, in Tik, the Hamiltonian H is represented by a so called pseudo-Jacobi ma¬ 
trix, that is, a Jacobi matrix pre or pos multiplied by the involution matrix 
diag(l, -1,1, -1,...), 


/3k -|- |fc| -|-1 

--f^y\k\ + 1 

0 

0 

7^1^! + 1 

jSk -|- |/c| -|- 3 

-7V2(|fc|+2) 

0 

0 

7V2(|fc| + 2) 

f3k -\-\k\-\-5 

-7^3(1^!+ 3) 

0 

0 

-f^y3{\k\ + 3) 

13k J- |/c| -|- 7 ... 






For simplicity, this matrix has been denoted by the same symbol H. It is clear 
that H is unitarily similar to its transpose 

j^T _ gi(7r/2)Ao^g-i(7r/2)no 
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The diagonalization of this matrix provides an alternative procedure to the 
previous approach. 

For completeness, it may be in order to observer that it is also possible to 
represent the generator ab by 


"o 

\/\k\ + 1 

0 

0 

0 

0 

^2(1*1 + 2) 

0 

0 

0 

0 

\/3{\k\ + 3) ... 

0 

0 

0 

0 






the generator a*b* by 


A'_ 


and the generator —a*a — bb* by 


A' 


/ T 


1 

1 

1 

0 

0 

0 

0 

-\k\-3 

0 

0 

0 

0 

1 

1 

0 

• 0 

0 

0 

-\k\-7 ... 






The commutation relations [A'_^ A'j^] = Aq, [Aq, ^4^ = 2A'_^, [Aq, A'_] = —2A'_^ 
which characterize the sm(1, 1) algebra, are clearly satished. However, in this case, 
the role of the lowest weight state is played by a highest weight state, such that 
= 0 . 

Let V be the linear space spanned by Aq, and let 


Hq = H - I3kl = Aq + 7(H+ - H_). 


Obviously, Hq ^ V. Let 


[Hq,V] := {HoX-xHo-.xeV}. 

Consider the linear operator [Hq, V] : V ^ V. We look for xA^ +yA_ + zAq such 
that 


[Hq, {xAj^ + yA_ + zAq)] — \{xAj^ + yA_ + zAq). 

The eigenvector xA+ + yA_ + zAq is easily determined by the secular equation 


2 0 -27 
0 -2 -27 
—7 —7 0 

The eigenvalues and respective eigenvectors of this 3x3 matrix are given by 

T 



X 


X 


y 

= A 

y 


z 


z 


2 VTT 


T 


i + yrr^ 


7 


7 ’ 1 + ^1 + 72’ 


11 
























-2/rT7. + 


7 


-i + yrr^’ 


and, of course, 0, ( 7 , — 7 ,1)^. Let us consider now the matrices 


B+ = 

B_ = 
Bo = 


7 f 1+ + 72 


2 vTT^ 

7 

2 vTT^ 

1 


7 


A+ + 


7 


i + ^TT^ 


A- — An 


-i + ^TTy 

7 


A+ + 


7 


-1 +vTT^ 


y4_ + Aq 


(Ao + 7(^+ “ ^-)) • 


VTT^ 

The following CR’s of the sm(1, 1) algebra are satished 

= Ro, [-Bo,-B+] = 2R+, = —2-B-, 

suggesting that 

= a/ 1 + 7^{|fc| + 1, |fc| + 3, |/c| + 5,...}, 


( 10 ) 


in agreement with the previous results in Section 01 We remark that Bq 7 ^ 
Rg, -B+ 7 ^ B*_. The operator (-B_) is said to be a raising (lowering) opera¬ 
tor. That is, if T is an eigenvector of Rg, so that RgT = AT, then R+T is an 
eigenvector associated with an upward shifted eigenvalue, SgR+T = (A-|-2)i?+T. 
Similarly, i?_T 7 ^ 0 is an eigenvector associated with a downward shifted eigen¬ 
value, 5 g 5 _T = (A - 2)i?_T. Next, we observe that 



satishes i?_Tg = 0, as may be easily checked. It follows that the spectrum of Bq 
is bounded from below. A laborious but straightforward computation shows that 
the Casimir operator reduces to 

C ■=Bl- 2{B+B_ + B_B+) = Bl- 2Bo - AB+B. = {k‘^ - 1)1. 


Thus, 

implying that 
so that 

This conhrms (fTOj) . 


CTg = (5g2 - 2Ro)To = (fc2 - l)To, 

-BofE'o = (1^1 + 1)'I^0! 
o'(-Bo) = {|^| + 1) 1^1 + 3, \k\ -|- 5,...}. 
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6 Discussion 


In Section 2, the diagonalization of non-Hermitian Hamiltonians with real 
eigenvalues, which are expressed as quadratic combinations of bosonic operators, 
is briefly analyzed. It is shown that, quite generally, Hamiltonians of this class 
are diagonalizable in terms of dynamical pseudo-bosons, which are determined 
by the EMM. In Section 3, a non self-adjoint Hamiltonian which is expressed 
as a linear combination of sm( 1, 1) generators, is investigated. Its eigenvalues 
and eigenvectors have been determined with the help of a matrix T of size 4 
that is real Hermitian and is determined by the EMM. The investigated Hamil¬ 
tonian has a complete system of eigenvectors expressed in terms of the creation 
and annihilation operators of pseudo-bosons, and is orthogonal to the system of 
eigenvectors of the adjoint Hamiltonian. Complete vector systems constructed 
in terms of boson creation operators acting on the associated vacuum state are 
obtained. Inhnite matrix representations of the Hamiltonian in such systems are 
presented. 

It would be interesting to consider more general Hamiltonians of the investi¬ 
gated type. A challenging problem would be to analyze the existence of inhnite 
dimensional versions in the spirit of Theorem 12.11 
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